Extraction of coherent structures in a rotating turbulent flow experiment 
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The discrete wavelet transform (DWT) and discrete wavelet packet transform (DWPT) are used 
to extract and study the dynamics of coherent structures in a turbulent rotating fluid. Three- 
dimensional (3D) turbulence is generated by strong pumping through tubes at the bottom of a 
rotating tank (48.4 cm high, 39.4 cm diameter). This flow evolves toward two-dimensional (2D) 
turbulence with increasing height in the tank. Particle Image Velocimetry (PIV) measurements on 
the quasi-2D flow reveal many long-lived coherent vortices with a wide range of sizes. The vorticity 
field exhibits vortex birth, merger, scattering, and destruction. We separate the flow into a low- 
entropy "coherent" and a high-entropy "incoherent" component by thresholding the coefficients of 
the DWT and DWPT of the vorticity field. Similar thresholdings using the Fourier transform and 
JPEG compression together with the Okubo- Weiss criterion are also tested for comparison. We 
find that the DWT and DWPT yield similar results and are much more efficient at representing 
the total flow than a Fourier-based method. Only about 3% of the large-amplitude coefficients of 
the DWT and DWPT are necessary to represent the coherent component and preserve the vorticity 
probability distribution function (PDF), transport properties, and spatial and temporal correlations. 
The remaining small amplitude coefficients represent the incoherent component, which has near 
Gaussian vorticity PDF, contains no coherent structures, rapidly loses correlation in time, and does 
not contribute significantly to the transport properties of the flow. This suggests that one can 
describe and simulate such turbulent flow using a relatively small number of wavelet or wavelet 
packet modes. 

PACS numbers: PACS numbers: 47.27.-i 47.32.-y 47.32.Cc 



I. INTRODUCTION 

Large-scale ordered coherent motions occur in a wide 
variety of turbulent flows despite existing in a rapidly 
fluctuating background turbulence 0, . These "coher- 
ent structures" , which are associated with localized re- 
gions of concentrated vorticity, persist for times that are 
long compared to an eddy turnover time. The impor- 
tance of coherent structures in turbulence has become 
recognized through the use of flow visualization. Exam- 
ples of coherent structures identified in turbulent flows 
include hairpin vortices in boundary layer turbulence Q , 
plumes in turbulent convection, and vortices in turbu- 
lent shear flows. Coherent structures play a major role 
in the transport of mass and momentum, thus affecting 
transport, drag, and dissipation in turbulent flows. Due 
to their long lifetimes, coherent structures in the atmo- 
sphere strongly influence the exchange of heat, moisture 
and nutrients between different locations. In industry, 
the prediction and control of transport, drag, and tur- 
bulence is important in many processes Q , and in some 
cases flows can be modified through the control of coher- 
ent structures [EE1- 

The emergence of coherent vortices is especially strik- 
ing in geostrophic turbulence, where the Coriolis force 
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plays a dominant role 0, II]- Under the effect of rota- 
tion and stratification, geophysical flows develop large 
and robust coherent structures which can be identified 
and tracked for time much longer than their characteris- 
tic turnover time. Examples of such flows are high and 
low pressure systems, large vortical structures that are 
formed in Gulf Stream meander, Mediterranean eddies 
(Meddies), and the Earth's jet stream; all can be ob- 
served and tracked by satellite imaging of the atmosphere 
or the ocean surface (e.g. see satellite imaging web- 
sites E3 ) ■ Large coherent structures are not limited 
to the Earth; the atmospheres of other planets also re- 
veal structures such as Neptune's dark spot and Jupiter's 
zones, belts and Great Red Spot. 

Three-dimensional (3D) turbulence subjected to strong 
rotation develops columnar vortical structures aligned 
with the rotation axis, as observed in the present experi- 
ment and in previous laboratory experiments [ill Il2l IT3 , 
IbH ] . Correlation in this direction becomes large when the 
Coriolis force becomes large compared to inertial forces, 
and the flow proceeds towards a quasi-2D state. This 
two-dimensionalization allows energy to proceed toward 
larger scales through the inverse energy cascade and to- 
ward smaller scales through the forward enstrophy cas- 
cade , as observed in simulation 0, 0, and ex- 
periments on rotating |l2l [13J and non-rotating quasi-2D 
flows 19]. The cascades of energy and enstrophy lead 
to a spontaneous appearance of intense localized coher- 
ent vortices containing most of the enstrophy of the flow. 
Vorticity filaments outside of those coherent structures 
are distorted and advected by the velocity field induced 
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by the vortices. Thus, a rotating turbulent fluid flow 
can organize itself into large-scale coherent structures 
that are often long-lived compared to dissipative time 
scales. The presence of long-lived coherent structures in 
turbulent rotatin g flo w has been observed in both exper- 
iment 0, [H milll and in simulation [HI El]. These 
structures are larger than the scale of the forcing. The 
large-scale coherent structures break the homogeneity of 
the flow and are thought to dominate the flow dynam- 
ics. The long lifetimes and spatial extent allow coherent 
structures to pla y a significant active role in transport 
processes |20ll2ll |. 

Due to the dynamical importance of coherent struc- 
tures, any analysis of flow containing these structures 
should take into account their existence. One approach 
to analyze the coherent structures and the dynamics of 
such flow is to partition it into regions with different dy- 
namical properties. Okubo |22j ] derived a criterion to 
separate flow into a region where strain dominates (hy- 
perbolic region) and a region where vorticity dominates 
(elliptic region). The same criterion was later re-derived 
by Weiss [23] and is now known as the Okubo- Weiss cri- 
terion. The criterion has been widely used to analyze 
numerical simulations of 2D turbulence. However, as 
pointed out by Basdevant and Philipovitch [24|, the va- 
lidity of the criterion's key assumption is restricted to 
the core of the vortices that correspond to the strongest 
elliptic regions. This limitation reduces the applicability 
of a decomposition using this criterion. 

Another method that has been found useful in an- 
alyzing flow fields is proper orthogonal decomposition 
(POD). This projects a field onto a set of orthonormal ba- 
sis functions where successive eigenvectors are obtained 
by numerically maximizing the amount of energy corre- 
sponding to that eigenmode [25] . [POD is known by other 
names, including Karhunen-Loeve decomposition, princi- 
pal components analysis (PCA), and singular value de- 
composition (SVD); the basis functions are also known as 
empirical eigenfunctions and empirical orthogonal func- 
tions [23, p. 86.] Linear combinations of the basis func- 
tions ideally correspond to the coherent structures in the 
flow. However, if the ensemble of fields is homogeneous, 
the basis functions become Fourier modes [25|. Further- 
more, a single mode of the POD corresponds to a full 
field with structures in a particular spatial arrangement. 
If the structures in the field are not stationary, as in 
the case of our present flow, many modes will therefore 
be necessary to track the different spatial configurations. 
Indeed, application of POD to our data resulted in the 
extraction of the large scale, low amplitude mean flow 
structure of our field, rather than the continually moving 
intense coherent structures. We therefore do not include 
POD in this paper. 

This paper separates a flow into coherent and incoher- 
ent components using wavelet transforms j2(| to extract 
localized features at different spatial scales. The most im- 
portant advantage of the wavelet representation over the 
more usual Fourier representation is the localization of 



the basis functions. A Fourier analysis is not well suited 
to pick out localized features such as intense vortices. 
The basis functions of a Fourier transform are localized 
in wavenumber space and hence spread out over the en- 
tire domain in physical space. The basis functions of 
the wavelet transform consist of dilates and translates of 
a "mother" wavelet, which contains multiple frequencies 
and has compact support (non-zero values only inside a 
finite interval) in physical space. The basis functions are 
well localized in both physical and wavenumber space; 
hence only a few coefficients suffice to describe localized 
features of a signal. The coefficients of the wavelet trans- 
form contain not only amplitude but also scale and po- 
sition of the basis elements. Thus the coefficients can be 
used to track the size and location of features that are 
well correlated with the wavelet bases (e.g. p?! l^l. 

The discrete wavelet transform (DWT) [2(j has been 
found to be well suited to analyze intermittent signals 
and systems containing localized features such as the in- 
tense vortices that occur in turbulence IHIlS Ip. Farge 
has extensively applied wavelets to the ana lysi s and com- 
putation of turbulent flows [H H El |H EH H| . 

The discrete wavelet packet transform (DWPT) [35| is 
a generalization of the DWT; the possible wavelet packet 
basis elements are a larger set which include spatial mod- 
ulation of the mother wavelet. The advantage of the 
DWPT is that the choice of basis is adaptable to the 
signal to be analyzed. 

Turbulent flow can be considered as a superposition 
of large-scale coherent motions, "fine-scale" incoherent 
turbulence, and a mean flow with interaction between 
the three constituents 1361 . In numerical simulations of 
2D turbulence [23, IH l3l| . and more recently 3D turbu- 
lence |33| . the coherent and incoherent turbulent back- 
ground components have been separated using wavelet 
based decompositions operating on the vorticity field. 
The coherent part, represented by only a small fraction of 
the coefficients, retained the total flow dynamics and sta- 
tistical properties, while the incoherent part represented 
no significant contribution to the flow properties. This 
separation of the flow into two dynamically different com- 
ponents suggests that the computational complexity of 
turbulent flows could be reduced in simulations with co- 
herent structures interacting with a statistically modelled 
incoherent background |3l| . The application, however, 
has been heretofore primarily limited to results obtained 
from numerical simulations. 

In this paper, we use the wavelet technique to analyze 
our Particle Image Velocimetry (PIV) data on a rapidly 
rotating turbulent flow. Section [H] describes the exper- 
imental system. Section ITU describes the resulting flow 
fields obtained in the experiment. Section llVI presents 
the techniques used to decompose the vorticity field into 
coherent and incoherent components. Section Ivl presents 
the results obtained by applying the method to measure- 
ments on rotating turbulent flow. The conclusions are 
discussed in Section IVTI 



3 



II. EXPERIMENT 
A. Instrumentation and measurements 

An acrylic cylinder (48.4 cm tall, 39.4 cm inner diame- 
ter) is fit inside a square acrylic tank (40 x 40 cm cross sec- 
tion, 60 cm tall) that has a transparent lid (Fig. 37] . 
The tank is filled with distilled water at 24 ± 1°C (p = 
0.998 g/cm 3 , v = 9.5 xl0~ 3 cm 2 /s). The tank and the 
data acquisition computer are mounted on a table that 
can be rotated up to 1.0 Hz. 
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FIG. 1: (a) A schematic of the apparatus. The height of 
the laser sheet is adjustable. Shown in the inset is a close-up 
a single tee. (b) Horizontal cross section of the distribution 
of pumping sources (open circles) and sinks (black dots) (po- 
sitions are to scale; see text for relative size of sources and 
sinks and tank diameter scale), (c) Overhead image of tank 
bottom, showing the tees. Overlayed are lines indicating the 
orientation of the tees. 

Fluid is injected at the bottom of the tank by pumping 
water through a hexagonal array of tubes with grid spac- 
ing 2.3 cm; there are 192 sources and 61 sinks, as shown 
in Fig. The source tubes are tee-shaped with a 2.5 
cm shaft and 2.1 cm horizontal top (0.079 cm inside di- 
ameter) [see inset of Fig. da)]. Each tee is screwed into 
the distributor. This design allows us to easily change 
the type of connector or change the forcing geometry by 
blocking or modifying sources or sinks. The tee source 
geometry was chosen to produce horizontal velocities and 



to minimize the vertical motions directly driven by the 
forcing. Further, the orientation of the tees was chosen 
to minimize cooperation between neighboring tees; how- 
ever, it was not possible to eliminate completely a mean 
flow. The sinks are 1.27 cm diameter holes in the bot- 
tom of the distributor [Fig. db)]- The forcing system 
is versatile and allows us to inject nearly homogeneous 
turbulence at a scale much smaller than the system size. 

Flow rates can range up to 1400 cm 3 /s, correspond- 
ing to flow velocities up to 8 m/s out of the tees. This 
corresponds to a Reynolds number of order 10 5 based on 
the jet velocity and the grid spacing of the tees. The 
Rossby number, U/2f2L, based upon the same scales is 
order 50 near the bottom of the tank, near the forcing. 
However, the turbulence decays away from the forcing, 
with increasing height, and becomes more 2D due to the 
influence of rotation. 

The data presented are for 0.4 Hz rotation rate and 
flow rates of 426 cm 3 /s, producing 2.4 m/s velocity jets 
at the forcing tees. This corresponds to Reynolds and 
Rossby numbers of about 6 x 10 4 and 20 respectively, 
near the forcing. Near the top the RMS flow velocity 
is about 2 cm/s and the maximal characteristic velocity 
length scale (twice the e-folding length of the velocity 
spatial correlation function) is about 5 cm; hence near 
the top the Reynolds number is about 1000. The Rossby 
number calculated from characteristic scales is 0.08, while 
the local averaged Rossby number calculated from the 
RMS vorticity, (wjjms/2^), at such a flow condition is 
0.3. The depth of our system thus allows us to observe 
more 3D (Rossby « 20) or more 2D (Rossby « 0.3) flow 
without changing control parameters, such as rotation or 
forcing. 

The water is seeded with polystyrene (p — 1.067 
g/cm 3 ) spherical particles with diameters in the range 
90-106 /im. The small mismatch in density results in 
a sedimentation terminal velocity of 4 x 10 -2 cm/s in 
the absence of flow, which is insignificant compared to 
the measured flow velocities. The density mismatch also 
causes a lag in the response of the particles in regions of 
large acceleration (33, H3 • We estimate the largest lag to 
be a 1% difference between the particle velocity and the 
flow velocity in the vortices at the top of the tank. 

A 395 mW, 673 nm diode laser with attached light- 
sheet optics (from Lasiris) illuminates the particles in 
horizontal planes for flow visualization and PIV measure- 
ments. It is also possible to rotate the laser and optics by 
90° to illuminate vertical planes for measurements of the 
vertical velocities. The thickness of the sheet is about 1 
cm and varies less than 10% across the diameter of the 
tank. The laser is fixed to a carriage which allows us 
to adjust the vertical position of the light-sheet without 
changing the pumping rate or system rotation rate. 

Particles are imaged in the rotating frame with a 
CCD camera (1004 x 1004 pixels, 30 frames per second) 
mounted above the tank. The images are grabbed and 
stored into a memory buffer by a computer on the ro- 
tating table. A second, lower resolution camera with an 
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analog output allows us to view the flow when not grab- 
bing digital images. 

The laser pulses were timed so that pairs of pulses were 
imaged in successive frames of the CCD camera, which 
had a dead time of 120 /is between consecutive frames. 
The pulse duration was typically 10 ms with a 2 ms dark 
interval between pulses. 

The analog signal from the second camera and power 
for the equipment on the table are sent through slip rings. 
The entire experiment is controlled remotely with a wire- 
less ethernet connection and remote access software (PC 
Anywhere) by a computer in the non-rotating laboratory 
frame. This allows us to take multiple data sets, adjust 
the pumping rate, and adjust the position of the light 
sheet without having to stop the table. 

To determine the two-dimensional projection of the 
velocity field we use a type of Digital Particle Image 
Velocimetry (DPIV) known as Correlation Image Ve- 
locimetry (CIV) | 40ll4l| . The CIV algorithm uses cross- 
correlations over small interrogation regions between a 
pair of consecutive images to find particle displacements 
over known time intervals to obtain velocities. The CIV 
algorithm also allows for possible rotations and defor- 
mations of the interrogation region to take into account 
nonlinear particle motions in between images. The re- 
sulting velocity field is interpolated to a uniform regular 
grid by a cubic spline interpolation. The vorticity fields 
are calculated from the coefficients of the spline fit. 

We can achieve a maximum data rate of 15 velocity 
fields per second (using two images to calculate each ve- 
locity field). This time resolution is sufficient for our 
flow conditions. To capture the dynamics of the flow, 
data were taken for 20 s periods (limited by memory 
buffer size), resulting in 300 fields in a single sequence. 
To capture statistical information, data were taken at 30 
s intervals (longer than measured decay time for our flow) 
over (in most cases) runs of 50 minutes in duration. This 
was sufficient to achieve convergence of statistical quan- 
tities. The resulting velocity and vorticity fields are on a 
128x128 grid with a spatial resolution of 0.30 cm when 
the whole tank is imaged and 0.19 cm when we zoom-in. 
Testing the CIV algorithm against test particle images 
in known flow fields 0) C3 > we found that the velocity 
fields determined by the algorithm have about 2% RMS 
uncertainty. 



B. Decay time 

Characteristic decay times for our flow were measured 
by three different types of experiments. In each experi- 
ment, the measured decay time was taken as the e-folding 
relaxation time of the mean kinetic energy in the flow. 

(i) Laminar (no forcing) spin-down experiments were 
conducted with the tees removed and replaced with a 
flat horizontal boundary. The tank was subjected to a 
sudden 10% decrease in rotation rate after the flow had 
reached solid-body rotation (no motion in the rotating 



tank frame). The predicted decay time of motions in 
laminar rotating flows with rigid, flat horizontal bound- 
aries subject to small step changes in rotation is given by 
the Ekman dissipation time, r = H / (2 V z/fi) , where H is 
the depth of the fluid, v is the kinematic viscosity and 
f2 = 2nf is the angular frequency of the container |44| . 
The corresponding decay time for the energy in the flow 
is t = H/(4^/ vfl). For our closed cylindrical tank with- 
out topography of depth H = 48.4 cm, v = 0.095 cm 2 /s, 
and a rotation rate of 0.4 Hz, this gives a decay time of 
the energy as 78 s; the measured time was 64 s. 

(ii) Laminar spin-down experiments were conducted in 
the same way as in (i) , except that the tees were installed 
in the bottom of the tank. The characteristic decay time 
measured with the tees installed in the bottom of the 
tank for our flow was 18±3 s. The reduced laminar spin 
down time is due to the tees, which cause extra drag 
and secondary circulations that quickly bring the fluid 
to solid-body rotation. 

(iii) Turbulent decay experiments were conducted by 
abruptly shutting off the forcing. For these experiments, 
the flow was allowed to reach a steady turbulent state 
(typically many decay times) under a constant pump- 
ing rate of 426 cm 3 /s and rotation rate of 0.4 Hz before 
abruptly turning off the forcing. The measured decay 
time of 13±3 s is long compared to the typical vortex 
turnover time of 1 s. 



C. Passive scalar advection 

To study the transport and mixing properties of our 
flow we examine numerically the motion of passive scalar 
point particles and passive scalar fields in the velocity 
fields obtained from the experiment. The velocity fields 
that we measure are a two-dimensional projection of a 
three-dimensional incompressible flow field. Therefore, 
they have a non-zero divergence and do not satisfy any 
fluid dynamical equation of motion. Nonetheless, they 
can give us useful information regarding the transport 
properties of the flow. 

Initial locations are chosen for the point particles, the 
positions of the particles are updated corresponding to 
the velocity fields such that x„ + i = u„Ai + x n where x ra 
is the position of the particle at time-step n, and u„ is the 
velocity of the flow field at time-step n at the location of 
the particle. Each velocity field is interpolated in space 
by a cubic-spline to calculate the field at the location 
of a given particle. The timestep, At, is chosen by the 
Courant condition |45j . which avoids particles jumping 
over grid points or going too far at a given iteration. For 
our data this condition means At < 0.02 s < 1/15 s; 
therefore, fields must also be interpolated in time. The 
measured velocity fields vary slowly in time compared 
to our temporal resolution 1/15 s (measured correlation 
time is ~ 2.4 s; see section lV Afl . We therefore justify the 
use of a cubic spline interpolation in time to achieve At 
< 1/15 s (typically At = 0.001 s). 
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We also examine the time evolution of a passive scalar 
field advected in the experimental velocity field by nu- 
merically integrating the advection-diffusion equation, 



TT7 = -U-Vc+KflV'c, 

at 



(1) 



where c = c(x, t) is the passive scalar concentration 
and kd is a diffusion coefficient, chosen as necessary for 
numerical stability of the solution. The values of 
used correspond to Schmidt numbers (kd/V) near 0.05. 
The numerical integration is performed using a pseudo- 
spectral method in polar coordinates based upon meth- 
ods given in [46j. The grid is 128 Chebyshev modes in the 
radial direction and 256 Fourier modes in the azimuthal 
direction. The velocity field obtained by the CIV mea- 
surement is interpolated via a cubic spline onto the simu- 
lation grid. The numbers of radial and azimuthal modes 
were chosen so as to not under-resolve our velocity fields 
on the nonuniform polar grid. 

The singularity at the origin was avoided by choosing 
Chebyshev modes for the radial direction. The colloca- 
tion points for Chebyshev modes have unit spacing on 
a circle, Xi = cos(wi/(N — f)), where Xi is the i th grid- 
point, and N the total number of gridpoints. The grid- 
points are clustered at the boundaries and sparse in the 
center of the domain. By defining the radial coordinate 
r 6 [—1,1] and using only the r > values, the radial 
grid points are sparse near r = and clustered at the 
boundary r = R. For an even number of modes the ori- 
gin, r = 0, is skipped. However, the azimuthal grid is 
still dense at the origin, which requires us to use a very 
short time step At w 2 x f 0~ 4 s to avoid numerical insta- 
bility. The experimentally determined velocity fields are 
then interpolated in time by the method described above 
for the tracer particle simulation. 

The diffusive term is calculated implicitly by a Crank- 
Nicholson scheme, separately for the radial and az- 
imuthal directions. The advection term is calculated ex- 
plicitly by a predictor-corrector scheme using a third- 
order Adams-Bashforth step followed by a fourth-order 
Adams-Moulton step [45| . 

To perform the numerical analysis on the decomposed 
fields we construct a velocity field from a vorticity field. 
We make the assumption of 2D flow that the total vor- 
ticity is given by the measured vertical vorticity, w = lo z . 
We then use the 2D streamfunction-vorticity relation 
V 2 * = —cu z and solve Poisson's equation for the stream- 
function \1/ using Matlab's PDE solver. The derivatives 
of the streamfunction are then used to calculate the ve- 
locity field by the relation V x (^z) = u, where z is the 
unit vector in the vertical direction. 

The assumption above is clearly not valid close to the 
source tees, where the flow is 3D. However, close to the 
top of the tank the assumption becomes valid as the flow 
is quasi-2D. To test the reconstruction and the validity 
of the 2D approximation, we compare the original mea- 
sured velocity fields to velocity fields reconstructed from 
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FIG. 2: (color on-line) Close-up of a vortex in our tank near 
the boundary (bold line). The original field is represented 
by black vectors and the reconstructed field by light-colored 
vectors. The original and reconstructed fields are indistin- 
guishable in the regions of uniform flow. The largest vectors 
correspond to 6.7 cm/s. The center of the tank is at (x,y) = 
(19.2 cm, 18.9 cm). 



the vorticity fields. A region of an original velocity field 
and the reconstructed field is shown in Fig.|2 The RMS 
difference in the magnitudes of the original and recon- 
structed fields is about 2%. The reconstruction calcula- 
tion does well in regions where there is a strong uniform 
flow. It does less well where there are large gradients in 
the velocity, in particular near vortices. However, on the 
whole the reconstructed velocity field follows the same 
behavior as the original field. 



III. VELOCITY AND VORTICITY FIELDS 
A. Transition to quasi-2D flow 

Near the forcing sources (tees) at the bottom of the 
tank, the flow is very turbulent (Reynolds number « 
6 x 10 4 ) and three dimensional (Rossby number 3> 1). 
However, moving vertically away and up the tank from 
the forcing, the turbulent velocities decay. The Reynolds 
number near the top of the tank where our data were 
collected is order 1000, based upon the RMS velocities 
of 2 cm/s and typical structure size of 5 cm (estimated 
from the velocity correlation function). The relative in- 
fluence of rotation also becomes larger so that the Rossby 
number becomes about 0.3. 

Fig. |21 compares the magnitude of the divergence and 
vertical vorticity of the flow in our tank at two different 
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FIG. 3: (color on-line) The PDF of the magnitude of vertical 
vorticity and divergence at two different heights in our tank, 
4 cm below the lid (solid curves) and 10.4 cm above the top 
of the tees (dashed curves). Inset: the ratio of the divergence 
PDF to the vorticity PDF, showing the increase in the relative 
magnitude of the divergence near the forcing. 

heights. The divergence in our flow fields is small rela- 
tive to the vorticity and 2S7 (ss 5 rad/s). The divergence 
field consists of small length-scales near the top and has 
a weak correlation with the vorticity field. The ratio of 
the RMS divergence to the RMS vorticity is 0.2, of the 
same order as the Rossby number. Near the bottom forc- 
ing, however, the divergence field becomes larger in am- 
plitude, length-scale, and more strongly correlates with 
structure in the vorticity field. The inset in Fig. f3] shows 
the ratio of the divergence to vorticity increases near the 
forcing, where the flow is more 3D and decreases near the 
top, where the flow is more 2D. 

The particles in a plane near the top of the tank are 
confined primarily to the plane, and the particle streaks 
follow persistent coherent jets, cyclones (a vortex with 
rotation in the same sense as the tank) , and anti-cyclones 
without crossing, similar to a 2D streamline flow [Fig. 0] 
(a) and (b)]. In contrast, the motions of the particles near 
the bottom of the tank are much less organized; particles 
pass through the plane frequently and the particle streaks 
can cross and exhibit small fluctuations due to increased 
three-dimensionality and turbulence. 

B. Persistent coherent structures 



flow [see Fig. 0|(d) for a close-up of an intense cyclone]. 
The anti-cyclones, however, have a typical core vorticity 
amplitude only a few times that of the background. 

The cyclones, anti-cyclones, and vorticity filaments are 
long-lived and are active dynamically. The vortices and 
vortex filaments can be tracked by eye as the flow evolves 
in time [Fig. 01(e)]. Vortices are shed from the wall and 
travel across the tank and interact with one another. Vor- 
tex filaments occasionally peel off of vortices, advect with 
the flow, and sometimes roll up to form new vortices. 
Structures may disappear by merging with other struc- 
tures or by stretching in a jet region between opposite 
sign vortices. Many of the small intense vortices with 
very short turn around times (< 1 s) live on the order 
of a characteristic decay time (~ 10 s) before shearing 
apart or merging with neighboring structures. 

Additionally, we see persistent structures such as the 
large cyclone (red) and anti-cyclone (blue) that appear 
in Fig. 0] Such structures often appear in a preferred 
location in the tank. Presumably such preferred loca- 
tions exist because of inhomogeneities in the forcing by 
the tees [cf. Fig. ^ ( c )]- A vortex can be kicked off its 
preferred location by a large perturbation when it inter- 
acts with neighboring vortices. If it begins to wander 
around the tank, it will generally disappear within a few 
characteristic decay times unless it returns to the pre- 
ferred location. If a structure moves off its preferred 
location and disappears, a new structure will typically 
form, replacing the pre-existing one within a few decay 
times. Some long-lived coherent structures occasionally 
persist throughout the entire duration of an experimen- 
tal run (~ 5000 vortex turnover times). Further, we ob- 
serve persistent structures as low as 5 cm above the tees. 
Observations closer to the tees are difficult because the 
three-dimensional turbulent flow rapidly moves particles 
into and out of the laser sheet. 

A wide range of spatial scales is visible in the veloc- 
ity and vorticity fields [Fig. 0] (c) and (d)] . The largest 
features, approximately 10 cm in size, are coherent vor- 
tices that have a much larger amplitude than their sur- 
rounding region. There are also large amplitude vortex 
filaments that stretch up to 10 cm; these can be as thin 
as the grid resolution in the transverse direction. Small 
scale persistent structures less than 1 cm in size yet large 
in amplitude are also observed. The region between the 
various large amplitude coherent structures is occupied 
by relatively low amplitude vorticity [see the vorticity 
trace in Fig. 01(c)] . 



Localized coherent structures including large compact 
regions of intense vorticity and wispy filaments of intense 
vorticity are evident in Fig. 0] (c) and (d). We observe 
more cyclones (rotation in the same sense as the tank) 
than anti-cyclones, in accord with observations in previ- 
ous experiments on rotating turbulent flows [l2T.ll3|. The 
amplitudes of vorticity at the core of the intense cyclones 
are typically more than ten times that of the surrounding 



IV. DECOMPOSITION OF THE VORTICITY 
FIELD 

A. Wavelets and wavelet packets 

The discrete wavelet transform (DWT) is a multi- 
resolution analysis that successively decomposes the sig- 
nal into coefficients which encode coarse and fine details 
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FIG. 4: (color on-line) (a) [MOVIE Streak photo (400 ms exposure time) of the particle image fields in a horizontal plane 
4 cm below the lid of the tank, (b) Close-up of the boxed region (18.7 cm x 18.7 cm) showing a cyclone (circular closed particle 
streaks) and an anti-cyclone (elliptical closed particle streaks to the upper- left of the cyclone), (c) Vorticity field ui z with the 
value along the dashed line shown in the trace below the field. The lines indicating the flow structure would be streamlines if 
the velocity field were divergence free, which it is not. Arrows indicate flow direction, (d) Close-up of velocity and vorticity 
fields in a 8 cm x 8 cm region near the boundary. The longest velocity vector corresponds to 6.7 cm/s. The vorticity is 
indicated by the color map, where vorticity values are clipped at 5 s _1 to render visible weaker structures; the vorticity in the 
core of the strongest cyclone is 22 s _1 . (e) [MOVIE |48| l A sequence showing the time evolution of the vorticity field. 
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at successively lower resolution |35|. The basis elements 
of the transform %j) s p correspond to dilations and trans- 
lations of a mother wavelet function ip, where s is the 
scale (dilations) of the wavelet and p its position (trans- 
lations). Successive levels of the transform continue to 
split the coarse detail coefficients, effectively analyzing 
the signal at coarser and coarser resolution. 

The discrete wavelet packet transform (DWPT) is a 
generalization of the discrete wavelet transform. The ba- 
sis elements of the wavelet packet transform include, in 
addition to dilations and translations, spatial modulation 
of the mother wavelet at different resolutions. The basis 
elements ip s ,p,k take on an additional parameter k, which 
roughly corresponds to the modulation of the wavelet 
packet. In contrast to the DWT, the choice of basis of 
the DWPT is not unique j26j. The wavelet basis is con- 
tained within the possible choices of wavelet packet bases. 
To select the particular basis to use, a natural choice is 
the wavelet packet basis into which the coefficients of the 
transform most efficiently represent the signal. This is 
known as the "best basis". The best basis is typically 
calculated based upon the minimization of an effective 
entropy measure of the coefficients [351 , thus minimizing 
the "information cost" of the coefficients in the best ba- 
sis. The flexibility in basis choice of the DWPT allows the 
transform to adapt the basis to the particular signal be- 
ing analyzed. For example, if the signal contains regions 
of rapid fluctuations, the basis choice will reflect that by 
including more basis elements with high modulation. 

We use the DWT [O(N) operations] and DWPT 
[0(N \og 2 N) operations] on our experimentally obtained 
vorticity fields, as previous authors have done using 2D 
turbulent flow data from numerical simulations in a peri- 
odic square domain (e.g. [27], H3, OOI ) ■ We use the Matlab 
wavelet toolbox and the coiflet 12 (coif2 in Matlab no- 
tation) as the analyzing wavelet. The coiflet family of 
wavelets (see reference [26|. p. 258) has both compact 
support and can generate an orthogonal basis. These 
two properties allow one to select the localized features 
of the coherent structures and to treat the decomposition 
of the vorticity field as two orthogonal components. The 
choice of wavelet does not alter the results significantly 
as long as the basis is sufficiently smooth |3l| . 



B. Coherent structure extraction 

We use an algorithm for coherent structure extraction 
based upon a de- noising algorithm (e.g. chapter 11 of 
reference |35jh In the de- noising, the assumption is that 
the original signal can be represented by a few large am- 
plitude coefficients of an orthogonal transform using an 
appropriate set of basis functions, while the noise is con- 
tained in the many remaining coefficients of small am- 
plitude. The de-noising is then performed by applying 
a threshold to the resulting coefficients of the transform. 
Coefficients above the threshold amplitude are assumed 
to correspond to the signal while coefficients below the 



threshold correspond to the "noise" . The few large am- 
plitude coefficients are called the "coherent" coefficients, 
while the many small amplitude coefficients are "incoher- 
ent" . The coherent and incoherent parts of the vorticity 
are then reconstructed by the inverse transform. 

The choice of the threshold separating the coherent 
and incoherent parts of the signal is based upon a mea- 
sure of the number of significant coefficients N a , which is 
the theoretical dimension of the signal, defined by 

N (f) = e ff(/) (2) 

where H(f) is the entropy of a discrete signal / = [/,] 
(where [fi] is the set of discrete values of an arbitrary 
signal /). In our case, the fi's correspond to the dis- 
crete values of our vertical vorticity measurements. The 
entropy H(f) is defined as 

N 

H(f) = -J^PilogPi (3) 

i=i 

where pi = |/i| 2 /||/|| 2 is the normalized square modulus 
of the i th element of the signal, with N the number of 
elements and ||/|| 2 = J2i \fi\ 2 - N indicates how many 
of the largest coefficients should be retained to give an 
efficient, low entropy, representation of the signal. 

Thus, the decomposition algorithm is the following. 
We take the transform of an individual measured vortic- 
ity field. For the DWPT, the best basis is used. Then 
we find the number of significant coefficients N Q of the 
transformed vorticity field in the transform basis. The 
threshold is therefore based upon the value of the N^ h 
largest coefficient. Coefficients whose modulus is larger 
(smaller) than the threshold correspond to the coherent 
(incoherent) part of the vorticity field. We then take the 
inverse transform to get the coherent part in physical 
space. The incoherent field is reconstructed by subtract- 
ing the coherent from the total field. The process is then 
repeated for each of our vorticity fields. Our algorithm 
contains no adjustable parameters other than the selected 
wavelet and is based upon the assumption that the flow 
field has a low entropy component, corresponding to the 
coherent structures, and a high entropy component, cor- 
responding to an incoherent background. 

Our algorithm is different than the various iterative 
algorithms of Farge et al. [3l| and othe rs Wft. l35| . In 
the de-noising algorithm described in [27], l35j|. successive 
iterations extract a fraction of the coherent coefficients 
until a stopping criterion is met. The fraction of retained 
coefficients at each step and stopping criterion are ad- 
justable parameters. However, the basis of the DWPT 
can change at each iteration. The number of coefficients 
retained therefore loses its meaning since coefficients can 
be selected from completely different bases. Therefore, 
we only perform a single iteration without changing the 
basis for the DWPT. In the algorithm of Farge [31|, an 
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a priori assumption is made of Gaussian white noise in- 
coherent component superimposed upon a non-Gaussian 
component of coherent vortices. The threshold is then 
iteratively found to separate the two without adjustable 
parameters. 

C. FFT, JPEG compression, and Okubo- Weiss 
based techniques 

To examine the relative efficiency of the wavelet-based 
decompositions, we compare the DWPT and DWT to 
three other algorithms; one based upon a Fourier trans- 
form, a second on basic JPEG compression, and a third 
on the method developed by Okubo and Weiss |2^| • 

For the Fourier transform algorithm we threshold the 
coefficients of a 2-dimensional fast Fourier transform 
(FFT) of our vorticity. The coherent fields are then con- 
structed via an inverse FFT of the largest amplitude coef- 
ficients. The incoherent remainder fields are constructed 
by subtracting the coherent fields from the original fields. 

The basic JPEG algorithm consists of subdividing an 
image into 8x8 blocks and taking the discrete cosine 
transform over the sub-blocks. The coefficients are ar- 
ranged based upon the global average of the amplitudes 
among the sub-blocks. The average of the amplitudes is 
used to determine which modes to keep. After the thresh- 
old, each sub-block retains the same number of modes, 
each in the same position. Thus the number of modes 
retained must be a multiple of the number of sub-blocks 
which compose the image (for our f28 x 128 fields, this 
is 256). For more information on the JPEG image com- 
pression standard see reference |4(|, chapter 3. 

The Okubo- Weiss criterion splits the fields into ellip- 
tic and hyperbolic regions, dominated by vorticity and 
strain respectively. The regions dominated by vorticity 
and strain are then separated and taken to be the coher- 
ent and incoherent fields. For details on the Okubo- Weiss 
criterion see reference ^4J . 



V. RESULTS 
A. Decomposed vorticity fields 

The resulting coherent and incoherent vorticity fields 
obtained from the wavelet packet, wavelet, Fourier, and 
JPEG techniques are compared in Fig. |S| The Okubo- 
Weiss criterion was also tested for comparison with the 
other decomposition methods. The fields shown are con- 
structed with the same number of coefficients for the 
DWPT, DWT, and Fourier transforms, and the closest 
approximate number for JPEG. 

The coherent fields in Fig. appear to retain the large 
scale coherent structures [cf. Fig. 0Jc) and (e)]. The 
DWPT, DWT, and Fourier preserve the structure of the 
total vorticity field. The Okubo- Weiss criterion excises 
primarily the regions of large amplitude vorticity, which 



generally corresponds to the cores of vortices. However, it 
extracts only the cores and leaves behind the peripheries. 
The JPEG does poorly with so few coefficients; it barely 
picks out some of the stronger vortices in the field. 

For the same number of coefficients, 2.4%, the DWPT 
and DWT do a better job of extracting structure in 
the fields than the Fourier decomposition, which leaves 
more structure behind in the incoherent field (Fig. [5] 
[MOVIE H|]). The JPEG incoherent field is large ampli- 
tude because of the few coefficients retained in the coher- 
ent field. The incoherent fields resulting from the DWPT, 
DWT, and the Fourier are mostly devoid of large-scale 
structures and are much smaller in amplitude. The in- 
coherent fields are also poorly correlated in time and are 
devoid of large scale structures or any other feature which 
can be tracked by eye (Fig. [5] and associated movie). 

The Okubo- Weiss criterion selects the centers of the 
vortices where vorticity dominates strain. This accounts 
for a large portion of the enstrophy. It extracts about 
74% of the enstrophy in the vorticity dominated regions, 
which cover about 40% (not coefficients) of the flow. 
However, the method leaves behind holes in the vortic- 
ity field that act as coherent structures and are often 
surrounded by large values of vorticity. The resulting 
"incoherent" field thus retains coherent properties. The 
Okubo- Weiss criterion also does poorly in recovering the 
statistics of the total field, as shown in Table Q] Fur- 
thermore, the decomposition shown in Fig. [3] was per- 
formed without regard to the Okubo- Weiss validity cri- 
terion, which restricts the application of the criterion to 
the very centers of the vortices or regions of very large 
vorticity p4| . If the validity criterion is applied, much 
less is retained in the coherent field, and the incoherent 
field ends up with large values of vorticity. 

Statistical properties of the decompositions are shown 
in Tables HI and ITT1 Table [i] lists the results from setting 
the threshold of the various transforms based upon the 
number of significant coefficients calculated from equa- 
tion @. On average about 2-3% of the large-amplitude 
coefficients of the wavelet-based decompositions were re- 
tained in the coherent fields, which account for about 
85% of the total enstrophy of the flow. The Fourier re- 
tains roughly the same enstrophy but requires more co- 
efficients, about 4%. The JPEG decomposition retains 
only 67% of the total enstrophy with about 13% of the 
large-amplitude coefficients. 

Owing to the compact support of their basis functions, 
fewer coefficients were needed for the DWPT and DWT, 
based upon the entropy criterion in equation than 
for the Fourier and JPEG decompositions. Compression 
curves for enstrophy of the various decompositions ap- 
plied to our vorticity fields are shown in Fig. The 
DWPT and DWT both do similarly well at extracting 
the enstrophy in the fields with only a small number of 
coefficients. The DWPT probably does slightly better 
due to the adaptability of its basis. The Fourier based 
methods do not converge as rapidly as the wavelet-based 
methods. This is likely due to the non-locality of the ba- 
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DWPT DWT Fourier JPEG Okubo-Weiss 




FIG. 5: (color on-line) Decomposed vorticity fields, TOP: Coherent BOTTOM: Incoherent remainder fields, all for the same 
number of coefficients retained (except for Okubo-Weiss and approximate for JPEG). The colormap corresponds to the intensity 
of vorticity and is the same for all images and the same as used in Fig. H [MOVIE [H(|] 



Decomposition 




Wavelet packet 




Wavelet 




Fourier 




JPEG 


Okubo- 


-Weiss 


Quantity 


Total 


coherent-incoherent 


coherent-incoherent 


coherent-incoherent 


coherent- incoherent 


coherent-incohere: 


coefficients retained (%) 


100 


2.4 


97.6 


2.7 


97.3 


4.1 


95.9 


13.1 


86.9 


40.2 


59.8 


Enstrophy /field (s~ 2 X 10 4 ) 


4.33 


3.63 


0.63 


3.67 


0.62 


3.68 


0.56 


2.92 


1.33 


3.21 


1.12 


Vorticity skewness 


0.78 


0.80 


0.02 


0.85 


0.04 


0.62 


0.20 


0.53 


0.24 


1.22 


-0.07 


Vorticity kurtosis 


7.09 


7.03 


3.88 


7.58 


3.42 


5.37 


7.85 


5.27 


6.55 


11.81 


8.28 



TABLE I: Statistical properties of the decomposed vorticity fields using entropy criterion. 



sis functions when applied to a field which has localized 
structures and sharp features. 

The wavelet-based algorithms also do a better job pre- 
serving the skewness and kurtosis of the total vorticity 
in the coherent component, while the incoherent compo- 
nents are much closer to a Gaussian distribution (skew- 
ness = 0, kurtosis = 3). Table ILT1 displays the results 
from setting the threshold on the transforms so that each 
method retains 2.4% of the coefficients in the coherent 
component (since the coefficients JPEG must be a multi- 
ple of the number of 8 x 8 blocks, 256, this restriction is re- 
laxed). For the same number of retained coefficients, the 
DWT and DWPT clearly outperform the Fourier method 
in terms of preserving the statistics of the total vorticity 
field. 

Figure [7\ shows that the coherent fields from the 
wavelet packet and wavelet decompositions rapidly con- 
verge to the statistics of the total vorticity field. The co- 
herent fields preserve the non-Gaussianity of the vorticity 
PDF without having to extract many coefficients. The 
respective incoherent fields also rapidly converge to Gaus- 
sian statistics. This suggests that wavelets and wavelet 
packets do a good job of extracting the coherent struc- 
tures from the vorticity field. The near-Gaussianity of 
the incoherent field suggest that the transforms have left 



the remainder without coherent structure. 

The coherent components of the Fourier and JPEG de- 
compositions, however, converge to the statistics of the 
total fields much more slowly. To obtain the same level of 
fidelity, both the Fourier and JPEG must extract many 
more coefficients. This inefficiency would result in over- 
transformed fields. Further, the insets in Fig.[7|(b) show 
that the kurtosis of the incoherent component never con- 
verges to anything small, but increasingly deviates from 
Gaussian. This suggests that the Fourier decomposition 
is not separating the coherent structures from the back- 
ground. 

The probability distribution function (PDF) of vortic- 
ity has broad wings (Fig. which correspond to the 
large vorticity values that occur in the coherent vortices 
and vorticity filaments. The preference for cyclonic (posi- 
tive vorticity) structures over anti-cyclonic structures ap- 
pears as a large positive skewness (see also Table H]or ITU . 
While the PDF for the DWPT and DWT coherent vortic- 
ity field are nearly the same as that for the total vorticity 
field, the PDF for the corresponding incoherent field is 
narrower and more symmetric, indicating the lack of high 
intensity structures. 

Consider the enstrophy spectrum of the vorticity fields, 
as shown in Fig. [H] (b). It is clear from the enstrophy 
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Decomposition Wavelet packet Wavelet Fourier JPEG 

Quantity Total coherent-incoherent coherent-incoherent coherent-incoherent coherent-incoherent 



coefficients retained (%) 


100 


2.4 


97.6 


2.4 


97.6 


2.4 


97.6 


3.0 


97.0 


Enstrophy/field (s~ 2 X 10 4 ) 


4.33 


3.63 


0.63 


3.60 


0.69 


3.28 


0.92 


0.96 


3.18 


Vorticity skewness 


0.78 


0.80 


0.02 


0.86 


0.04 


0.51 


0.32 


0.06 


0.72 


Vorticity kurtosis 


7.09 


7.03 


3.88 


7.62 


3.47 


4.77 


7.85 


4.84 


7.08 



TABLE II: Statistical properties of the decomposed vorticity fields retaining the same number of modes, 2.4%, as the DWPT 
decomposition (approximate for JPEG). 



1.0 




0.2 



0.1 0.2 

Fraction of coefficients retained in coherent component 

FIG. 6: (color on-line) The percent enstrophy retained in the 
largest amplitude coefficients as a function of the number of 
coefficients kept. In order from best (most efficient) to worst: 
wavelet packet, wavelet, Fourier, JPEG. 



spectra of the total field that most of the enstrophy is 
contained in long wavelengths (small k) . This agrees well 
with the observation of the vorticity fields and the dom- 
inance of large structures with large amplitude vortic- 
ity. The enstrophy spectrum, shown in Fig. [SJ does not 
contain a well defined scaling region corresponding to a 
cascade of enstrophy. A cascade would not be expected 
for our quasi-2D flow, which is forced by the broad-band 
3D turbulence in the bottom of the tank rather than by 
injection of energy at a single well-defined wavenumber. 
The presence of large structures which extend well into 
the 3D region alone tells us that we must have a broad 
spectrum energy injection. 

The coherent field contains the same large structures 
that are present in the total field. Accordingly, the en- 
strophy spectrum of the coherent field matches that of 
the total field at long wavelength [Fig. |SJ (b)]. The inco- 
herent field contains negligible enstrophy at large scales, 
indicating the lack of large scale features. This is ap- 
parent in the vorticity decomposition (see Fig. [SJ) . The 
incoherent field has only small structures and is the dom- 
inant contribution to the enstrophy at large k. 

The wavelet-based and the Fourier decompositions also 
retain the spatial and temporal correlations in the coher- 



ent field, as Fig. El illustrates with the DWPT and DWT 
yielding almost the same results. The long time correla- 
tion is in part due to the presence of long lived coherent 
structures. In contrast, for the incoherent field the spa- 
tial and temporal correlations are short ranged, indicat- 
ing the absence of large scale and long-lived structures. 



B. Transport of passive scalar particles and 
concentration fields 

Paths computed for some passive scalar "particles" in 
the measured velocity fields are shown in Fig. ^| (a). 
Over the 20 s duration of the measurement, the particles 
are advected around a large portion of the tank. The 
motion of the particles depends on their location in the 
flow. For example, a particle can spend time caught in a 
vortex, following the motion of the vortex as it meanders 
slowly in a localized region of the tank [path in upper 
right of Fig. EH (a)] . Occasionally the particle will escape 
from a vortex and then be carried by a high velocity jet, 
which can transport a particle large distances in a short 
time. A particle may then be captured by a vortex [e.g. 
path ending in lower right of Fig. ^| (a)] . 

The resulting mean squared displacement of the tracer 
particles in the velocity fields obtained by experiment is 
shown in Fig. ^| (b) . It exhibits approximately t 2 scaling 
at short times, which looks "flight-like" 21]. At long 
times the finite size of the tank becomes significant and 
the scaling exponent becomes smaller. 

Both the wavelet and wavelet packet coherent fields 
disperse the tracer particles similarly around the tank; 
the scaling behavior of the mean squared displacement, 
not shown in this paper, is also the same. Differences 
in the paths are due to the sensitive dependence upon 
initial conditions of a turbulent flow field. 

We find a striking difference in the behavior of a nu- 
merically integrated passive scalar particle in the coher- 
ent and incoherent fields. The incoherent fields make 
no significant contribution to the transport properties; 
rather, particles are confined to a small region, as shown 
in the inset in Fig.^|(a). This is due to the rapid decor- 
relation in time of the incoherent fields. The rapid fluctu- 
ations cause the tracers to jiggle around and the scaling 
of the mean squared squared displacement of the trac- 
ers in the incoherent fields approaches that of a random 
walk. 
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The evolution of the passive scalar held in the mea- 
sured velocity held is similar to that for tracer particles. 
Results for the coherent and incoherent fields produced 
by the DWPT are shown in Fig. ^2 which illustrates 
the stretching and folding; the results for the total held 
are similar to that of the coherent held. By four sec- 
onds the scalar has already been signihcantly stretched 
by the velocity helds. There is no signihcant advection 
in the incoherent helds and is similar to results obtained 
with pure diffusion (no advection term) . This is expected 
from the short time and space correlations of the inco- 
herent helds (Fig. 0). Similar results were obtained for 
the DWT and Fourier methods. Our results agree well 
with those of Beta et al. j^il who, using the DWT, found 
that the coherent held is responsible for the mixing in a 
numerically simulated 2D turbulent flow. 

VI. DISCUSSION 

We find that flow which is strongly forced near the bot- 
tom of a deep rapidly rotating tank evolves with increas- 
ing height in the tank from 3D turbulence into rotation- 
dominated turbulence. This is due to the competition 
between the turbulent forcing which is 3D, and the ro- 
tation which tends to two-dimcnsionalizc the flow. Near 
the forcing the turbulent velocities are large and the flow 
is 3D. Away from the forcing the spatial decay of tur- 
bulence allows the rotation to become increasingly im- 
portant with increasing height. Motions of the particles 
in the flow are more organized in the horizontal plane 
near the top, consistent with two-dimensional columnar 
vortices. The how in the rotationally dominated (Rossby 
number ss 0.3) region of our tank has a strongly non- 
Gaussian vorticity PDF and exhibits properties of 2D 
turbulence, such as large-scale long-lived coherent vor- 
tices, vortex filamentation, and merger. 

The vertical component of vorticity describes the flow 
well near the top where the flow is quasi-2D as a con- 
sequence of the rotation. We are thus able to use the 
projected flow field as measured in a horizontal plane for 
analysis with wavelet and Fourier-based decompositions. 
We have shown that wavelet based transforms can be 
used to separate the coherent structure-containing quasi- 
2D turbulent flow into a non-Gaussian coherent compo- 
nent represented by as few as 3% of the large-amplitude, 
low entropy coefficients of the transform and a nearly 
Gaussian incoherent component, represented by the re- 
maining small-amplitude, high entropy coefficients. We 
are able to obtain this result without an a priori assump- 
tion of Gaussianity or non-Gaussianity of the two com- 
ponents. 

The discrete wavelet packet transform (DWPT) and 



the discrete wavelet transform (DWT) yield very similar 
results, despite the adaptability of the DWPT basis. The 
DWT is therefore made preferable by its faster compu- 
tation [O(N) versus 0(N \og 2 N) operations]. 

Because our flow helds contain compact structures, the 
localized basis functions of the DWPT and DWT out- 
perform Fourier and JPEG decompositions. Both the 
DWPT and DWT have more rapid convergence of the 
statistics of the extracted coherent component toward 
that of the total flow. The rapid convergence of the skew- 
ness and kurtosis of the vorticity PDF suggest that the 
wavelet based methods have efficiently captured the co- 
herent structures. The Fourier and JPEG converge much 
more slowly and are thus unable to efficiently capture the 
large higher moments of the vorticity PDF. This suggests 
that the Fourier and JPEG methods have not really ex- 
tracted the coherent structures despite the appealingly 
smooth visual appearance of the structures in the coher- 
ent held of the Fourier. Indeed, the superior performance 
of the DWT over the JPEG is exploited by the emerging 
next generation JPEG2000 image compression standard, 
which uses a bi-orthogonal DWT |49|, |52| . The incoher- 
ent remainder of both the DWPT and DWT converge 
rapidly towards Gaussian statistics, while incoherent re- 
mainders of the Fourier and JPEG do not converge to 
any value. 

The coherent components of the DWPT and DWT re- 
tain all of the properties of the total held, including the 
large-scale structures, shape of the vorticity PDF, long 
spatial and temporal correlations, and transport proper- 
ties. Further, the coherent component contains the large 
skewness and kurtosis of the PDF which are due to the 
coherent structures. In contrast, the incoherent remain- 
der has only small-scale short-lived features and does not 
contribute signihcantly to the transport. Thus in analy- 
sis of how dynamics and transport, it may be sufficient 
only to consider the coherent component. These results 
suggest that it is reasonable to reduce the computational 
complexity of turbulent flows by considering only the low- 
dimensional coherent structures, which interact with a 
statistically modelled incoherent background. 
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FIG. 7: (color on-line) Convergence of the skewness and kur- 
tosis of the vorticity PDFs for the various decompositions, 
plotted with quantities defined so that they approach unity 
as they converge. The vertical line indicates the fraction of 
coefficients retained in the wavelet packet decomposition, (a) 
Convergence of the coherent component. Plotted: kurtosis 
ratio: |fc co h/fetotal|, skewness ratio: [s co h/stotal|- (b) Conver- 
gence of the incoherent component. Plotted: kurtosis ratio: 
2 &mcoh/3, skewness ratio: 1 jsincoh/stotal 

I. Note that 

unity on the graph represents the respective statistic for a 
Gaussian distribution: skewness = 0, kurtosis = 3. The kur- 
tosis of the Fourier and JPEG incoherent field is large and 
off the main plot. The insets show the convergence behavior 
at large numbers of retained coefficients. The open circles on 
the curves correspond to the fraction of coefficients retained 
for the corresponding techniques using the entropy criteria as 
a threshold. 
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FIG. 8: (color on-line) (a) Vorticity PDFs of the total field, 
coherent field, and incoherent field, where the inset shows the 
incoherent PDFs plotted together with a Gaussian fit to the 
wavelet packet results, (b) Enstrophy spectra Z(k) of total, 
coherent, and incoherent fields, calculated on a square sub- 
section of our circular domain. The solid, dashed, and dash- 
dotted lines correspond to the wavelet packet, wavelet, and 
Fourier decompositions respectively. The k 1 line is the spec- 
tral slope for Gaussian white noise in 2D. The wavenumber is 
defined as k = 2n/L. 
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FIG. 9: (color on-line) (a) Space and (b) time correlations of the decomposed vorticity fields. The solid, dashed, and dash-dotted 
curves correspond to the DWPT, DWT, and Fourier decompositions. 
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FIG. 10: (color on-line) (a) Simulated paths of numerical 
tracer particles in the total field (solid curves) and in the 
coherent (dashed curves) and incoherent field (red) for the 
DWPT (total time 20 s). The inset shows an expanded view 
of the path in the incoherent field, (b) Mean squared displace- 
ments for the tracers in the total, coherent, and incoherent 
velocity fields. The exponents for the total and coherent fields 
are 1.97 while for the incoherent remainder the exponent is 
1.32 over the range indicated. Both results are for the wavelet 
packet. 
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FIG. 11: (color on-line) [MOVIE 0] Advection of a passive 
scalar field in the velocity field of the DWPT, as computed 
from the advection-diffusion equation Q. By eight seconds 
the scalar field has been significantly mixed by the coherent 
field, while it appears only to have been diffused in the inco- 
herent field. 



